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Summary 

Results obtained during the first year of an investigation of the 
structure of sprays are briefly described. Further details may be found 
in Refs. 1-3. 

The investigation involves both experimentation and analysis. Experi- 
mental objectives are to complete measurements of the structure of nonevapo- 
rating, evaporating and combusting .sprays for sufficiently well-defined 
boundary conditions to allow evaluation of models of these processes (Fig. 1), 
Analytical objectives are to begin model evaluation using both existing and 
the new data (Fig. 2). The results of the investigation have application to 
the development of rational design methods for aircraft combustion chambers 
and other devices involving spray combustion (Fig. 3). 

Major assumptions for the models are summarized in Fig. 4. The con- 
tinuous phase is treated using a k-e-g model of turbulence originally pro- 
posed by Lockwood and Naguib [4], which has been extensively calibrated for 
noncombusting and combusting single-phase flows during earlier work in this 
laboratory [ 5-7 ] . 

Three methods for treating the discrete phase are being considered: 

(1) a locally homogeneous flow (LHF) model, (2) a deterministic separated 
flow (DSF) model, and (3) a stochastic separated flow (SSF) model. The 
main properties of these models are summarized in Figs. 5-7. Infinitely 
fast interphase transport rates and local thermodynamic equilibrium are 
assumed for the LHF model (Fig. 5) — implying that both phases have the same 
temperature and velocity at each point in the flow. LHF models provide a 
useful limit for infinitely small particles or drops, but generally over- 
estimate the rate of development of practical sprays [5-7]. DSF models 
(Fig. 6) provide for finite interphase transport rates, but assume that 
interphase transport can be found by ignoring effects of turbulent fluctu- 
ations. Most spray models reported to date employ this approximation. The 
present SSF model (Fig. 7) adapts an approach originally proposed by Gosman 
and loannides [8]. In this case, particles or drops are assumed to interact 
with a succession of, turbulent eddies whose properties are determined by 
random sampling — given mean and fluctuating properties of the flow from the 
k-c-g model calculations. This involves computation of a statistically 
significant number of particle trajectories using Monte Carlo techniques. 
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Initial model evaluation employed data for dilute particle-laden jets — 
to avoid complications due to particle coalescence, particle collisions and 
polydisperse particle flows. Only a sample of the results is given here, 
cf. Refs. [1-3] for complete findings. 

The prescription for eddy properties used in the SSF model was cali- 
brated using theoretical results of Hinze [9] — similar to Gosman and 
loannides [8]. This analysis was for the dispersion of inf initely-small 
particles in a homogeneous and isotropic turbulent flow. The comparison 
between SSF predictions and the analytical result is illustrated in Fig, 8, 

Good agreement was achieved — fixing methods for estimating particle inter- 
actions with eddies. The SSF model predictions are compared with measure- 
ments of particle dispersion in a duct flow, reported by Snyder and Lumley 
[10], in Fig. 9 , The model is seen to provide encouraging predictions of 
effects of particle properties on rates of particle dispersion. 

The remaining comparisons between predictions and measurements consider 
particle-laden jets. LHF, DSF and SSF model predictions are compared with 
the measurements of Yuu et al, [11] in Fig. 10. The LHF and DSF models 
over- and under-estimate particle dispersion, while the SSF model is in good 
agreement with measurements. Furthermore, the DSF model indicates that the 
particles tend to concentrate near the centerline as axial distance increases — 
which is not observed. Comparison of the models with measurements of McComb 
and Salih [12,13], Laats and Frishman [14,15] and Levy and Lockwood [16] 
continues in Figs. 11-13. The SSF model yields satisfactory predictions for 
these flows, aside from possible effects of turbulence modulation and turbu- 
lence generation at high particle mass loadings [2,3]. This evaluation is 
not adequately definitive, however, due to uncertainties in initial conditions 
for the existing particle-laden jet data [1-3], 

Evaluation of the models is continuing using data from the present in- 
vestigation. A sketch of the test apparatus being used for noncombusting 
sprays is illustrated in Fig. 14. An air-atomizing injector sprays vertically 
downward along the centerline of a traversible screened enclosure. Experi- 
mental methods are summarized in Fig. 15. All techniques have been used in 
work to date, aside from. the LDA-visibility method for drop size and velocity 
measurements — which requires a different optical geometry and is being de- 
ferred until other measurements are complete. 

Experimental methods were established by satisfactory measurements of 
the properties of air jets, formed by the injector, with earlier work [5,6]. 
Tests were then conducted in two nonevaporating sprays having SMD of 87 and 
30 ym. LHF and SSF model predictions are compared with measurements of mean 
gas velocity and mean liquid flux, along the spray axis, in Figs. 16 and 17. 
There is no fundamental limitation in the use of the LHF model in dense 
regions of the spray; therefore, these predictions extend from the injector 
exit. The separated flow models, however, are limited to dilute regions of 
the spray; therefore, these calculations begin at x/d = 50 — where adequate 
initial conditions were available from the measurements. The prediction of 
the LHF model improves for the more finely atomized spray. Case 1, but is 
not very satisfactory. In contrast, the SSF model provides reasonably good 
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predictions of these measurements. Predicted and measured radial profiles 
of liquid mass flux are plotted in Fig. 18. Similar to the results for 
particle-laden jets, the DSF model yields excessive concentrations of the 
dispersed phase near the axis, since turbulent particle diffusion is ignored. 
The SSF model provides fair predictions of particle spread. However, in 
this case, the LHF model underestimates spread rates I This effect is well- 
known, involving enhanced spreading of particle-laden flows by turbulent 
diffusion for a certain range of particle inertial properties. The fact 
that the SSF model correctly predicts this trend is very encouraging. Pre- 
dicted and measured turbulence kinetic energy are illustrated in Fig. 19. 
While the LHF model underestimates the magnitude of k and the width of the 
flow, the SSF model provides satisfactory predictions. This, too, is 
encouraging, since predictions of k are an important element in estimating 
eddy properties for the SSF model. Good predictions of Reynolds stress 
were also obtained with the SSF model. Reynold stress depends on z pre- 
dictions — suggesting that this aspect of the eddy prescription is also 
adequate. 

The conclusions, to date, are summarized in Fig. 20. The main con- 
clusion is that the SSF model provides encouraging predictions for these 
multiphase flows — with minimal added empiricism. It will be most interesting 
to examine this methodology for evaporating and combusting sprays — where 
effects of concentration fluctuations must be considered along with velocity 
fluctuations. The main limitation of the evaluation, thus far, is adequate 
specification of initial conditions for the present data base. Tests during 
the next report period are designed to eliminate this deficiency (Fig. 21). 
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Objectives:’ Complete measurements of spray structure suitable 
for evaluation of models. 

Configuration: Axisymmetric spray or particle-laden jet in a 

quiescent environment. 

Specific Cases: i. Air jet (calibration). 

II. Particle-laden jet (sand particles), 
hi. Non-evaporating spray (low-volatility oil). 
IV. Evaporating spray (Freon-11). 

V. Combusting spray (h-Pentane). 


Figure 1.- Experimental Objectives. 


Objectives; Complete evaluation of typical models using both 
existing and new data. 

Models: i. Locally homogeneous flow (LHF)— Infinitely 

fast interphase transport rates. 

ii. Deterministic separated flow (DSF)— Finite 
interphase transport rates considering 
particle response to mean motion. 

III. Stochastic separated flow (SSF)— Finite 
interphase transport rates but particles 
respond to individual eddies. 

Evaluation: Gas jets^ solid-particle-laden jets, nonevaporating 
sprays, evaporating sprays and combusting sprays. 


Figure 2.- Analytical Objectives. 
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Potential contributions of the results are: 

1. The development of rational design methods for 

AIRCRAFT combustion CHAMBERS^ GAS-TURBINE 
COMBUSTORSy FURNACESy STRAT I F I ED-CHARGE I.C. 
ENGINES AND DiESEL ENGINES, 

2. Improved understanding of the properties of 

TURBULENT PART I CLE/DROP-LADEN FLOWS, 


Figure 3.- Applications of the investigation. 


1. Steady^ AXisvMMETRiCy large Reynolds NUMBERy boundary-layer flow, 

2. k-c-g TURBULENCE MODEL WHICH IS WELL-CALIBRATED FOR NONCOMBUSTING 
AND COMBUSTING SINGLE-PHASE JETS, 

3. Negligible kinetic energy and viscous dissipation of mean flow 
AND radiation. 

Exchange coefficients of all species and heat identical. 

5, Continuous-phase is in local thermodynamic equilibrium. 

6. DSF AND SSF ONLY; Dilute -particulate flow so effects of turbulence 

GENERATION AND DISSIPATION BY PARTICLESy PARTICLE COLLISIONSy AND 
ADJACENT-PARTICLE DISTURBANCES OF INTERPHASE TRANSPORT RATES ARE 
NEGLIGIBLE. 


Figure 4.- Major assumptions of the models. 
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Assumption: Infinitely fast interphase transport rates, i,e., 

PARTICLE AND CONTINUOUS PHASE VELOCITIES AND 
temperatures are IDENTICAL AND LOCAL THERMODYNAMIC 
EQUILIBRIUM INCLUDES BOTH PHASES. 



Notes: i. 

II. 

in. 


Correct limit for infinitely small particles. 

Maximum particle response to turbulent fluctuations. 
Computation equivalent to single-phase flow. 


Figure 5.- Properties of the LHF model. 


Assumption: Finite interphase transport rates with particles 
responding to mean motion. 



pi 


pi 


^ U, 


ip 


T 

^G. 


Notes: i. Particle dispersion ignored— only valid for 

"large" particles. 

II. Typical approach in current spray models. 

III. Eulerian calculation for continuous phase with 

DISTRIBUTED SOURCE TERMS FROM PARTICLE INTERACTIONS. 

IV. LaGRANGIAN calculation of PARTICLE TRAJECTORIES. 

Figure 6.- Properties of the DSF model. 
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Assumption: Finite interphase transport rates with particles 

INTERACTING WITH A SUCCESSION OF INDIVIDUAL EDDIES 
WHOSE PROPERTIES ARE FOUND BY RANDOM SAMPLING OF 
LOCAL TURBULENCE PROPERTIES. 



, (2k/3)’''^) 

t|pUF(?, g’ 
= G.jpDF (7 


Notes: i. Maximum particle displacement and time of eddy 

• p 3/4 i3/2y . _ I 

INTERACTION ARE: *-e ' ^ ' Lg/(2k/3) 

II. Provides predictions of fluctuating particle 

PROPERTIES AND TURBULENT PARTICLE DISPERSION. 


III. Computation similar to DSF model--Monte Carlo 

TECHNIQUE TO FIND PARTICLE TRAJECTORIES. 

Figure 7.- Properties of the SSF model. 



Figure 8.- SSF predictions for the dis- 
persion of infinitely small particles 
in a homogeneous isotropic flow (ana- 
lytical results from Hinze [9]). 


20 






0 100 200 300 400 

t(ms) 


Figure 9.-. Predicted and measured particle dispersion in a 
uniform grid-generated turbulent flow (measurements from 
Snyder and Lumley [lO]). 



Figure 10.- Comparison of LHF, DSF, and SSF 
predictions of particle dispersion with 
the measurements of Yuu, et al. [ll]. 
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r/x 

Figure 13.- Predicted and measured mean 
and fluctuating particle velocities 
(data from Levy and Lockwood [16]). 



Figure 14.- Sketch of the experimental 
apparatus. 
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Measurement 


Technique 


Mean and fluctuating gas and 

PARTICLE VELOCITIES 

Drop (particle)— size distributions 

Drop-size distributions 

Drop size and velocity 
Mean liquid (particle) flux 

Mean temperature 
Mean composition 


Laser-Doppler anemometer (LDA) 

Fraunhofer diffraction (monitor- 
ing nonevaporating flows only) 

Slide impaction 

LDA-visibility method 

Isokinetic sampling and filtering 
(nonevaporating flows only) 

Shielded, fine-wire thermocouple 

Isokinetic sampling and analysis 
with gas chromatograph 


Figure 15.- Summary of experimental methods. 



x/d 


Figure 16.- Predicted and measured mean gas 
velocities along the axis of nonevaporating 
sprays . 
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Figure 17. Predicted and measured mean liquid 
mass flux along the axis of nonevaporating 
sprays . 











1. Model Performance: 

I. LHF AND DSF MODELS ARE VALID AT LIMITS OF SMALL AND LARGE 
PARTICLES^ BUT DATA BASE (AND PROBABLY SPRAYS IN GENERAL) 
INCLUDED FEW RESULTS AT THESE CONDITIONS. 

II. SSF MODEL (PARTICULARLY THE NSSF VERSION) YIELDED ENCOURAGING 
RESULTS— WITH MINIMAL EMPIRICISM. 

2. Decisive model evaluation was not achieved due to inadequate 

SPECIFICATION OF INITIAL CONDITIONS FOR MOST OF THE DATA BASE. 

3. Effects attributable to turbulence generation and modulation 

WERE OBSERVED^ HOWEVER^ THE EFFECTS WERE SMALL SINCE. DATA BASE 
WAS LIMITED TO DILUTE SPRAYS (VOID FRACTION > 99%). 


Figure 20.- Conclusions of the investigation* 


Task 

Status 

Model Development 

Completed for noncombusting 
flows. Combusting flows next. 

Measurements 


Single-Phase Jets 

Completed. 

Particle-Laden Jets 

Completed, aside from mean and 
fluctuating particle velocities. 

Nonevaporating Sprays 

Completed^ aside from drop size 
and velocity measurements. 

Evaporating Sprays 

Completed^ aside from mean and 
fluctuating gas velocities (in 
progress) and drop size and 
VELOCITY measurements. 

Combusting Sprays 

Next. 


Figure 21.- Status of the investigation. 
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